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The Tunka Radio Extension (Tunka-Rex) is an array of 44 radio antenna stations, distributed over 
3km^, constituting a radio detector for air showers with an energy threshold around lO'^eV. It 
is an extension to Tunka-133, an air-Cherenkov detector in Siberia, which is used as an exter¬ 
nal trigger for Tunka-Rex and provides a reliable reconstruction of energy and shower maximum. 
Each antenna station consists of two perpendicularly aligned active antennas, called SALLAs. An 
antenna calibration of the SALLA with a commercial reference source enables us to reconstruct 
the detected radio signal on an absolute scale. Since the same reference source was used for the 
calibration of LOPES and, in a calibration campaign in 2014, also for LOEAR, these three exper¬ 
iments now have a consistent calibration and, therefore, absolute scale. This was a key ingredient 
to resolve a longer standing contradiction between measurements of two calibrated experiments. 

We will present how the calibration was performed and compare radio measurements of air show¬ 
ers from Tunka-Rex to model calculations with the radio simulation code CoREAS, conhrming it 
within the scale uncertainty of the calibration of 18%. 
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1. Introduction 

To overcome restrictions of commonly used detectors for cosmic ray air showers, new detector 
types and detection principles are being investigated since the inception of air shower detection. 
One promising candidate for a next generation air shower detector, or part of a hybrid detector, is 
the radio detection technique. It is based on the radio emission of air showers: Mainly due to the 
geomagnetic deflection of charged particles, a short electromagnetic pulse in the radio frequency 
range is emitted [1]. A small, but still significant contribution to the pulse from the Askary an 
effect [2], originating from a varying net charge, was experimentally confirmed lately [3, 4, 5]. 

The resulting pulses can be detected with an antenna array. The main advantages of this 
method are its almost full duty cycle and the indicated high sensitivity to air shower parameters. 
Because of the high transparency of the atmosphere for radio frequencies, it is also possibly well 
suited for high observation depths or highly inclined showers. 

The first generation of digital radio arrays [6, 7] proved their general feasibility and sensitivity 
to air-shower parameters, while also fueling progress in the theoretical understanding of the radio 
emission. Contemporary radio arrays now have to either show their competitiveness to established 
detection principles or refine the methods and technology to reach it. 

The Tunka Radio Extension (Tunka-Rex) is a digital radio-antenna array in the Tunka valley 
in Siberia, Russia. It is co-located with Tunka-133 [8], a photomultiplier array for the detection 
of cosmic-ray air showers via their Cherenkov emission in air. The detector concept of Tunka- 
Rex is to record the radio signal of a shower in coincidence with the air-Cherenkov signal, by 
triggering the antenna array with Tunka-133. Thereby obtaining a reliable trigger, as well as the 
shower parameters from an independent measurement with known precision to compare the radio 
reconstruction to. However, Tunka-133 is limited in its operation time to moonless nights during 
winter. To overcome this limitation, in near future, Tunka-Rex will also be triggered by Tunka- 
Grande [9], a recently installed scintillator-based detector for cosmic-ray air showers at the site, 
capable of operation around-the-clock. 

One goal of Tunka-Rex is to obtain a calibrated measurement of the radio signal. This enables 
further testing and improvement of the modeling of the radio signal, as well as determining its ab¬ 
solute scale for comparison to other radio experiments. For this purpose we performed an antenna 
calibration with a reference source. 

In this paper we describe the Tunka-Rex antenna station and its calibration. Finally we com¬ 
pare measured radio amplitudes from air-shower events to simulated ones with the simulation soft¬ 
ware CoREAS. 

2. Detector setup 

The layout of Tunka-Rex is depicted in Fig.l. Its central part consists of 19 antenna stations, 
distributed over 1 km^. These antennas are currently triggered by Tunka-133. There is one antenna 
connected to each of the 19 clusters of Tunka-133. A cluster consists of a group of 7 PMTs in a 
hexagonal pattern with common housing for electronics and data acquisition (DAQ) in its center. 
Together with 6 more distant antenna stations Tunka-Rex covers 3 km^, but these outer antennas are 
only relevant for the highest energies, since the antenna spacing decreases from roughly 200 m to 
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Figure 1: The layout of Tunka-133 and Tunka-Rex in its different stages. Tunka-133 is organized in hexag¬ 
onal clusters. Each cluster is equipped with one Tunka-Rex antenna station and a scintillator station with 
another antenna station connected to it. 

500 m in this outer region. With the Tunka-Grande scintillator extension each cluster was extended 
by a scintillator station with another antenna station, totaling 44 antennas, for the whole array. 
Once Tunka-Grande is fully comissioned, it will trigger the whole array. 

Each Tunka antenna station consists of two channels, each with a SALLA [10] antenna, 
aligned perpendicularly to each other. By reconstructing the incoming direction of the radio signal 
from the relative arrival time, the two channels enable the reconstruction of the electric field vector. 

Each SAEEA is directly connected to a low-noise-amplifier (ENA), in the upper antenna box 
(see Eig. 2), enhancing the signal by about 22 dB. Erom there the signal is transferred via coaxial 
cable to the cluster center in 20 m distance. The cluster center is heated to ensure stable operation 
for the contained electronics in cold winter nights. There, the signal is filtered to the frequency 
range 35-76 MHz and enhanced by another 32 dB by a filter-amplifier. Einally, fhe so far analog 
signal is digitized af 200Ms/s wifh 12 bif depfh, and a « 5 ps frace confaining 1024 samples is 
recorded around fhe frigger time of Tunka-133. 

3. Antenna calibration 

The defermine fhe elecfric field on an absolufe scale from fhe ADC counfs of fhe frace, fhe 
impacf of fhe full hardware chain has fo be understood and inverted. Erom a formal poinf of view, 
fhe impacf of a linear sysfem on an incoming signal is described by fhe convolution of fhe incoming 
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Figure 2: Setup of the Tunka-Rex antenna calibration in Karlsruhe. A crane holds a reference source above 
the antenna aligned from the ground using strings. With the known emitted field strength of the reference 
source, the characteristics of the SALLA can be determined. 

signal with its response function. For an antenna, this response function depends on the incoming 
direction. For dipole like antennas, the vector effective length H is commonly used as response 
function and especially convenient if the polarization of the signal is of interest. In case of an 
antenna the incoming signal is an electric field E and the outgoing signal is a voltage V. In the 
frequency domain the convolution becomes a simple point-wise multiplication 

V{v)=Hiv)-E{v). (3.1) 

Provided two measurements of V from the two perpendicular antenna channels with linearly inde¬ 
pendent sensitivity to the polarization, this equation can be inverted to reconstruct all three compo¬ 
nents of E, if the incoming direction and H is known. While the incoming direction in each event 
is determined from the arrival times, assuming the incoming direction of the shower and the radio 
signal coincide, H remains the same for all events and has to be determined only once. 

To obtain H, we performed a calibration of the full analog hardware chain, including the 
antenna, with a reference source. Therefore, a complete Tunka-Rex antenna station was deployed 
at the Karlsruhe Institute of Technology, Germany. On the Tunka site, the station is plugged into 
the existing DAQ, which was not available in Karlsruhe, so we had to rely on a oscilloscope for 
data acquisition. This makes the calibration almost end-to-end, except for the DAQ, containing the 
analog-to-digital converter, which was calibrated independently. 

The reference source we used was the same as the one used for the calibration of LOPES [11] 
and recently LOFAR [12]. The available calibration of the reference source was performed under 
free-held conditions, with a high impact from ground rehections. Therefore, we obtained new cali¬ 
bration values for the reference source for free-space conditions, which are also used for an update 
of the LOPES results [13]. Due to the use of the same reference source, these three experiments 
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now have the same amplitude scale and can be compared neglecting the scale uncertainty due to 
the calibration of the reference source. 

The reference source is the commercially available biconical antenna and signal generator 
VSQ1000+DPA4000 by Schaffner Electrotest GmbH (now Teseq). It produces a signal linearly 
polarized along the antenna. In the time domain the signal is a comb with one puls per ps. In the 
frequency domain this results in a comb, consisting of peaks with 1 MHz spacing. The absolute 
height of the peaks, i.e. the amplitudes at r^/ = 10 m for the respective frequencies are given by 
the manufacturer with an accuracy of 16% (one sigma), making up for the dominant uncertainty in 
the amplitude reconstruction. 

For the calibration, the source is placed at different zenith angles above the antenna in >10m 
distance to the antenna by a crane, standing in 20 m distance, holding the source with a wooden 
spacer. With a differential GPS the distance r between the reference source and the SALLA are 
determined with a precision of several 10 cm. The source is aligned via strings from ground along 
the polarization axis of maximum sensitivity of one of the SALLAs. So the one channel receives 
maximal power, while the other one receives minimal power. We can obtain H{v) via: 


Hiv) 


vjy) 

E{v) 


r 

- P- 

rref 


(3.2) 


E{v) is the absolute value of the electric field strength of the reference source. To obtain the 
oolarization axis of maximum sensitivity p, an antenna simulation of the passive antenna structure 
is performed with the NEC2 code. It turns out, that the sensitivity of the SAEEA closely follows 
the pattern of a horizontal dipole in the antenna arc plane. 

Measuring the response of the analog hardware separately with a network analyzer [14], we 
can subtract it from the total response to extract the VEE of the bare antenna. In Fig. 3, the simu¬ 
lated VEE is shown against the measured one. The relative deviation of measured from simulated 
VEE over different zenith angles is 12%, averaged over all relevant frequencies, used as an esti¬ 
mated for the precision of the directivity model. The scale deviation is 6%, which is in the expected 
order due to the calibration uncertainty of the reference source of 16%. 

For the calibration we normalize the simulated antenna pattern frequency-wise to the mean of 
all calibration measurements at different zenith angles. The total scale uncertainty is 18%, resulting 
from the squared sums of 16% from the uncertainty of the calibration of the reference source, 6% 
from the temperature during the calibration and 4% from positioning and alignment of the reference 
source. The uncertainty on the absolute amplitude reconstruction on antenna level is 22%. Besides 
the 18% on the scale, it includes 12% from the antenna model and some smaller contributions from 
temperature, antenna alignment and uncertainties in the production and deployment of the antenna 
stations (details in [15]). 


4. Amplitude comparison 

After the calibration we compare model calculations of the radio emission to the actual mea¬ 
sured amplitudes. Therefore, we take the Tunka-133 data set of the first season, 2012/2013, where 
the reconstructed primary energy is above lO'^ eV and an zenith angle below 50°. 
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Figure 3: Measured vector effective lengths of the bare SALLA (markers) and simulated ones (lines). The 
left-right seperation is only for illustration, since all curves are mirror-symmetric. 


We then cut the frequency range to the realized bandwidth of 35-76 MHz and remove some 
single frequency entries from the spectrum at every 5 MHz, where we frequently observe interfer¬ 
ences. Then the full hardware response is inverted. A noise level is determined by the RMS of 
a 500 ns window before the signal window and a signal S by the peak of the Hilbert envelope in 
the signal window. As radio event only those qualify which have at least 3 stations with a signal- 
to-noise ratio (S^/N^) of at least 10, which enables the reconstruction of the shower axis from the 
relative arrival times of the pulses. The shower axis then has to agree with the on reconstructed by 
Tunka-133 within 5°. 

This leaves 92 events in the data set for which we perform simulations of the radio signal with 
CoREAS [16]. The events are simulated for proton and iron primaries with initial parameters from 
the Tunka-133 reconstruction. Since the depth of the shower maximum cannot be directly set in the 
simulation, but emerges from other parameters and shower fluctuations, we repeat the simulations 
several times and take only simulations into account, which happen to have their shower maximum 
closer than 30g/cm^ to the Tunka-133 reconstruction. 

The remaining events undergo a full detector simulation and afterwards the same reconstruc¬ 
tion algorithm and cuts are used as for the real data. 20 iron and 28 proton events pass all re¬ 
quirements and are used for comparing measurements and simulations. Since every event contains 
multiple stations, there are 83 and 124 amplitude measurements for iron and proton, respectively. 

In Fig. 4 these ampltiudes are shown versus the measured ones. The scale for iron and proton 
are off by a factor of 0.96 and 0.84 respectively, within the scale uncertainty of the calibration of 
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Figure 4: Amplitudes simulated with CoREAS versus measurements from actual showers. For the simula¬ 
tion initial parameters from Tunka-133 where used. Error bars indicate the energy ucertainty of Tunka-133 
for the simulation and the uncertainty from background for the measurement. 

18%. Event-to-event and shower-to-shower uncertainties average out for the mean over several 
antennas from different events. 

5. Summary 

Tunka-Rex is a digital radio array for the detection of air showers in Siberia. Besides deter¬ 
mining the precision and overall benefit of a radio detector in combination with an air-Cherenkov 
detector, its goal is to determine the absolute scale of radio emission. Therefore, we performed 
an absolute calibration with a commercially available reference source, reaching an accuracy for 
an individual amplitude measurement of 22%, including 18% uncertainty on the absolute scale. 
The same reference source was also used by LOPES and LOEAR. These three experiments now 
have a consistent amplitude calibration. Comparing the measured amplitudes to the predictions of 
the simulation code CoREAS, we find agreement within the uncertainties. Therefore, Tunka-Rex 
confirms the applicability of CoREAS for the simulation of absolute amplitudes of radio emission 
by air showers. It also provides a measure of the absolute scale of radio amplitudes in its frequency 
band. 
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